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Abstract 

We address, by means of ab-initio calculations, the origin of the correlation that has been ob- 
served experimentally between the chemisorption energy of CO on nanoscale Cu(001) supported 
films and quantum-size effects. The calculated chemisorption energy shows systematic oscillations, 
as a function of film thickness, with a periodicity corresponding to that of quantum-well states at 
f crossing the Fermi energy. We explain this trend based on the oscillations, with film thickness, 
of the decay length on the vacuum side of the quantum-well states at the Fermi energy. Contrary 
to previous suggestions, we find that the actual oscillations with film thickness of the density of 
states per atom of the film at the Fermi energy cannot account for the observed trend in the 
chemisorption energy. 

PACS numbers: 73.20.At, 73.21.Fg, 82.65. +r, 68.43.-h 
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I. INTRODUCTION 



Understanding and controlling the structure-size dependence of the chemical activity in 
materials with reduced dimensionality is a longstanding target of surface science, whose 
practical significance stems from potential applications in catalysis, gas sensing, and anti- 
corrosion. A recent exciting development in this area is the observation, in metal films a 
few atomic layers thick, of a correlation between quantum-size effects and properties such 
as chemisorption^ and surface oxidation.-"- 

Danese et al.,— in particular, have observed experimentally systematic oscillations in the 
desorption temperature of CO adsorbates on ultrathin epitaxial Cu(OOl) films on fcc-Fe(OOl), 
as a function of film thickness. Such Fe-supported films give rise experimentally to well 
resolved series of quantum-well states (QWS's).- 1 ^ The largest CO desorption temperatures 
were found to coincide with Cu film thicknesses at which a QWS crosses the Fermi energy, 
Ep, in inverse photo-emission (IPE) spectra taken at normal incidence.- 

Typically, modifications in chemisorption properties of metal surfaces are achieved by 
monolayer or submonolayer deposition of a second metal, forming a surface alloy, and ex- 
ploiting related local surface mechanisms such as site blocking or charge exchanges Using 
thin-film QWS's instead is an interesting new approach to tailor surface chemisorption prop- 
erties. In spite of its interest, however, the precise mechanism and key parameter behind the 
QWS-related modifications of the chemisorption properties are not yet fully understood. 

Very recently, quantum-size effects have also been reported in the initial oxidation rate 
of ultrathin Mg, Al, and Pb supported films.-"- The effect was initially suggested to be due 
to periodic oscillations, with film thickness, in the magnitude of the density of states (DOS) 
of the films at the Fermi energy.-^^ The DOS per atom of the film is a parameter which 
is often invoked in model reactivity theories, 10 and which may account for the observed 
reactivity trends of the Pb films^ — although other parameters have also been invoked.— 1 ^ 
In the case of the Mg films, however, it was observed later on that the actual variation in the 
DOS per atom of the film at Ep could not simply account for the order-of-magnitude change 
in the initial oxidation rate of the films.— Moreover, for the Al films, significant differences 
were observed between the DOS and reactivity trends as a function of film thickness.- 

In the cases of the Mg and Al films, for which the precursor adsorption mode of the 
O2 molecule on the surface should be physisorption, the key parameter responsible for the 
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changes in the initial oxidation rate was then proposed to be the decay length in vacuum, 
A, of the electronic local density of states of the film at the Fermi energy.—^ Modifications 
in A can be expected to have a direct exponential influence on the electron transfer rate by 
resonant tunneling, which is believed to control the initial sticking of the O2 molecules on 
such surfaces.-^ The same argument, however, based on tunneling processes, does not hold 
in the case of CO on Cu(OOl), where the molecule is clearly chemisorbed and the modulated 
value is the chemisorption energy. In this case, the key factor responsible for the observed 
changes is still an open issue. 

In this work, we investigate by means of ab initio calculations, the correlation between 
the chemisorption energy of CO molecules on Cu(OOl) films and quantum-size effects. The 
observed systematic oscillations in the chemisorption energy can be understood in terms of 
periodic modifications, with film thickness, in the decay length of the QWS's at the Fermi 
level. The paper is organized as follows. In Section [Til we briefly present the calculation 
method and parameters used. In Section IIII Al we examine the electronic spectra of free- 
stranding and fcc-Fe-supported Cu(OOl) films with thicknesses in the range 3-12 monolayers 
(ML). In Section [III Bl we examine the CO chemisorption energy on the corresponding films. 
The trends and microscopic origin of the correlation between electron quantum-size effects 
and chemisorption energies are analyzed in Section nil CI Our conclusions are summarized 
in Section ITVl 

II. COMPUTATIONAL DETAILS 

The calculations were carried out within the generalized gradient approximation (GGA) 
to density-functional theory (DFT), using the Perdew-Burke-Ernzerhof (PBE) exchange- 
correlation functional.— We employed the pseudopotential-plane-wave method, as imple- 
mented in the PWscf code of the QUANTUM-ESPRESSO distribution. 17 For carbon, oxy- 
gen, iron, and copper, we used the Rabe-Rappe-Kaxiras-Joannopoulos ultrasoft pseudopo- 
tentials of the PWscf library. The semicore Cu 3d states were treated as valence states. The 
nonlinear core correction 7 to the exchange-correlation potential was used for Cu and Fe. 

In the experimental study by Danese et al.,— the Cu films were grown on a closely latticed 
matched template formed by 5-ML of fee Fe(OOl) deposited on a Cu(OOl) substrate. In our 
study, we considered both free-standing and fcc-Fe-supported Cu(OOl) films, with/without 
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the CO adsorbates. All calculations involving the Cu films on Fe were spin-polarized calcu- 
lations. The systems were modeled using slab geometries in supercells. The length of the 
supercell was set to 85 A in all cases. A kinetic-energy cutoff of 30 Ry was used for the 
plane-wave expansion of the Kohn-Sham orbitals. 

For the supported Cu films, we considered slabs containing 5 Fe(001) ML plus 3 to 12 
Cu(001) ML deposited on one side. The Cu theoretical lattice constant of 3.68 A was 
used to construct the slabs (the experimental value is 3.61 A); we neglected the effect of 
the small lattice mismatch between fcc-Fe and Cu.— To evaluate the work functions of the 
Cu(001) supported films, we employed larger, symmetric slabs of n-ML Cu/5-ML Fe/n-ML 
Cu (3 < n < 12). The smallest vacuum region used was 31 A, corresponding to a 29-ML slab. 
For the CO covered films, we used asymmetric slabs with CO adsorbed on one side of the 
slab. We considered the cases of the Cu(001)-(1 x l)-CO and Cu(001)-c(2 x 2)-CO surfaces, 
corresponding to CO coverages = 1 ML and O = 0.5 ML, respectively. Experimentally, 
the desorption measurements in Ref. Qwere performed for the Cu(001)-c(2 x 2)-CO surface. 

The adsorption site of the CO molecule is known from experiment to be the on-top 
site.— This is in agreement with our PBE-GGA calculations, and consistent with the re- 
sults of previous GGA calculations using the same type of pseudopotentials.— We note, 
however, that although our GGA calculations yield the correct adsorption site for CO on 
Cu(OOl), this is not so in general within GGA for other transition metal surfaces - where 
the underestimation of the CO gap gives rise to a different site preference with respect to 
experiment.— 1 ^ We therefore also checked that changing the CO gap (with a GGA + U 



scheme, and Uco — 0.75 eV, as in Ref. |22|) does not affect (within 0.1 meV) the calculated 
variations of the chemisorption energy with film thickness. 24 

The CO molecule is vertical on Cu(001), the experimental Cu-CO bonding distance is 
d(Cu-C) = 1.92 A and the C-0 distance is d(C-0)= 1.13 A. 20 From calculations for the 
Cu(001)-c(2 x 2)-CO surface, using a 12-ML thick Cu film, we find d(Cu-C) = 1.9 A and 
d(C-0)= 1.15 A, in good agreement with the experimental values 2 ^ and with previous GGA 
results for thinner Cu films. 21 We kept these distances fixed in the remaining part of our 
study. We also elected to use the same frozen Cu (Cu/Fe) slab geometry for the CO covered 
and uncovered Cu films, as our main purpose here is to understand the correlation between 
the CO-Cu(OOl) chemisorption energy and quantum-size effects, which are of electronic 
origin. 
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We evaluated the CO chemisorption energy as: 

E ch (CO) = —(E s i ab+C o — E s i ab — Eco) (1) 

where E s i ab and E s i ab+C o are the total energies of the slab without and with the adsorbed 
carbon monoxyde, respectively, and Eco is the total energy of the CO molecule in the gas 
phase. For the calculation of the isolated carbon monoxide molecule, we have used a cubic 
cell with a parameter of a = 15 A. The corresponding equilibrium C-0 distance was found 
to be d(C-0)= 1.14 A. 

In the self-consistent calculations for the (1 x 1) [c(2 x 2)] surfaces, the integrations over 
the Brillouin zone were performed using a 24 x 24 x 1 [18 x 18 x 1] k-point grid centered 
at T. A Gaussian smearing of the electronic levels of 0.01 Ry was used to determine the 
Fermi energy. The relative values of the chemisorption energy were found to be converged 
to within ~ 1 meV with respect to the kinetic-energy cutoff, k-point grid, and vacuum size. 
For the calculations of the DOS of the film at the Fermi energy, we increased the k-point 
grid to (40 x 40 x 1), in order to ensure a numerical accuracy of 0.0004 eV" 1 on the DOS 
value per atom and spin. To calculate the local density of states at Ep, used to evaluate the 
decay length A, we employed a (48 x 48 x 1) k-point grid and a Gaussian smearing of 0.005 
Ry. We determined A from a fit, assuming an exponential decay of the local density of states 
at distances beyond ~ 2.4 A from the outermost atomic plane. The numerical accuracy on 
A was estimated as 0.002 A. 

III. RESULTS AND DISCUSSION 

A. Quantum-well states of the free-standing and supported Cu(OOl) films 

In FigJUa), we show the calculated energy levels of the QWS's with wavevector ku = 
for the free-standing Cu(001) films, as a function of film thickness. The unoccupied levels 
at F are the states probed in normal-incidence IPE experiments.- In Figfjjb), we also show 
the energy bands of the Cu bulk states with wavevectors perpendicular to the films, i.e., 
along the T— X line (A direction) of the bulk Brillouin zone. The QWS's of the films at f 
originate from these bulk states. 

The QWS's of interest near the Fermi level, in FigJT^a), with energies in the range [Ep— 1.5 
eV, E F + 1.5 eV], derive from the upper part of the bulk Cu 4sp band, in FigfjTJb), which 
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FIG. 1: (a) Calculated energies of the quantum-well states at f in the free-standing Cu(001) films 
as a function of film thickness, (b) Dispersion of the Cu bulk bands along the [001] direction (i.e., 
the T— X direction of the bulk Brillouin zone), perpendicular to the surface of the films. The zero 
of energy corresponds to the Fermi level. 

crosses the Fermi energy and has its maximum at the X point. The spectrum in the energy 
range [Ep — 5 eV, Ep — 1.5 eV] is dominated by the more localized Cu 3c/-band states. The 
lower part of the bulk Cu 4sp band (below Ep — 5 eV) produces another series of QWS's, 
in Figflja), located in the energy range [Ep — 9 eV, Ep — 5 eV]. 

In Fig. [21 we display the band structures of the 6-ML and 10-ML free-standing Cu(001) 
films along the high-symmetry lines of the surface Brillouin zone. The low- dispersion 3d 
bands can be easily recognized in the energy range [Ep — 5 eV, Ep — 1.5 eV]. We also note 
that the two series of quantum-well levels at F, originating from the top and bottom part of 
bulk Cu Asp band along the A direction, give rise to two characteristic series of parabolic-like 
subbands, in Fig. [21 near the surface Brillouin-zone center. 

In FigfjTa), one can observe that the energies, near Ep, of the QWS's of the film at f 
increase with increasing film thickness and cross the Fermi energy at 5 ML and at 11 ML. 
This is in good agreement with the IPE measurements, where the crossing occurs at 5 ML 
and at 10-11 ML.- The trend of increasing energy with increasing film thickness of the T 
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FIG. 2: Band structure of the 6-ML (left panel) and 10-ML (right panel) free-standing Cu(OOl) 
films along the high-symmetry lines of the surface Brillouin zone. The zero of energy corresponds 
to the Fermi level. 



QWS's (which is opposite to the trend of the QWS's near E F in the Mg(0001), Al(lll), and 
Pb(lll) films" ) results from the negative effective mass of the Asp band of bulk Cu, at the 
X point, in the A-direction. The top of the Cu Asp bulk band, at about 1.5 eV above Ep in 
FigfTJ^b), corresponds to the bottom of the "inverted" 1-D quantum well which determines 
the energies, near Ep, of the QWS's of the film at T. With increasing film thickness, the 
quantum-well levels become increasingly closer to the bottom of the inverted quantum well 
at ~ 1.5 eV above Ep. 

The periodicity of the crossing of Ep by T QWS's can be derived from the Bohr- 
Sommerfeld rule,— considering the parabolic-like part of the bulk Cu Asp band in FigfJJb) 
and its intersection with Ep. The empty section of the Cu Asp bulk band runs over 17 % 
of the Brillouin-zone X — T line, and therefore every 1/0.17 = 5.8 ML a quantum-well level 
of the film at T should cross the Fermi energy, according to the Bohr-Sommerfeld rule.— 
This is in good agreement with the crossing periodicity of 6 ML we find in Figfl^a), and 
consistent with the 5-6 ML crossing periodicity observed in the IPE spectra.- Our results 
are also in good agreement with previous DFT calculations for the T-state spectra of the 
free-standing Cu(001) films and their analysis.— 1 ^ 

As the experimental Cu films in Ref. |l| were grown on fee Fe(001), we also investigated 
the effect of the Fe substrate on the QWS spectra. In Fig. [31 we show the calculated partial 
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density of states at F of the Cu(OOl) films on the fcc-Fe(OOl) substrate, as a function of film 
thickness. 28 The partial density of states of the film was obtained by summing the spin- up 
and spin-down atomic-projected density of states at F of the Cu atoms in the film. The 
results in Fig. [3] show that the sequence of thicknesses at which the QWS peak maxima 
cross Ep remains exactly the same as for the free-standing films, namely 5 and 11 ML. 
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FIG. 3: Partial densities of states of the Cu(001) films at T, for the Fe-supported Cu films. The 
film thickness increases from 3 to 12 atomic layers (bottom to top curve). 



B. Quantum-size effects on the CO chemisorption energy 

The calculated CO chemisorption energies of the free-standing and Fe-supported Cu(OOl) 
films are shown in Fig. HJ for the two CO coverages = 1 ML and O = 0.5 ML. The 
behavior as a function of film thickness is rather similar for the two CO coverages. For the 
free-standing films (and for both = 1 ML and O = 0.5 ML), the chemisorption energy is 
largest at 3 ML, and displays a local minimum at 7 ML and a local maximum at 10 ML; 
for = 1 ML, in addition, a local minimum (maximum) is present at 4 (5) ML. In the case 
of the Fe-supported films, both O coverages give rise to the same trends: local maxima are 
found in the chemisorption energy at 4 and 10 ML, and a local minimum is present at 6 ML. 
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FIG. 4: Calculated chemisorption energy of the CO molecule for the free-standing (upper panels) 
and Fe-supported (lower panels) Cu(001) films, as a function of film thickness. Left panels show 
the results for the Cu(001)-(1 x l)-CO surface and right panels for the Cu(001)-c(2 x 2)-CO surface, 
corresponding to CO coverages of 1 ML and 0.5 ML, respectively. 



The presence of the Fe substrate modifies thus the trends of the chemisorption energy 
at low film thickness (below 8 ML), giving rise to a pronounced maximum at 4 ML and 
also shifting the local minimum from 7 ML (in the free-standing case) to 6 ML (for the 
supported films). At larger film thicknesses, instead, a local maximum is found at 10 ML 
in all cases. The presence of the substrate gives rise to a systematic oscillatory behavior 
in the chemisorption energy, which improves the agreement with the experimental trend. 1 
Experimentally, the local maxima in the chemisorption energy occur at 5 and 10-11 ML, 
and a local minimum is observed at 6 ML.i 

For the supported films, the oscillatory behavior of the chemisorption energy shows a 
correlation with the periodic crossing of Ep by a QWS at T. In fact, the theoretical maxima 
in the chemisorption [at 4 and 10 ML, in Fig. H](c) and (d)] are found to occur systematically 
1 ML before the crossing of Ep by a QWS (at 5 and 11 ML, in Fig. [2D- It may appear 
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surprising, at first sight, that the presence of the substrate restores a trend which one would 
expect to be intrinsically related to the properties of the film. It should not be forgotten, 
however, that effects other than those induced by the QWS near Ep can be expected to have 
a major influence on the chemisorption energy for free-standing films only a few monolayers 
thick. Indeed, the presence of a free Cu surface with broken bounds a few Cu ML away from 
the CO covered surface may be expected to strongly modify (enhance) the strength of the 
CO-Cu bound. 

We find that the chemisorption energy for G = 0.5 ML is considerably larger than that 
for = 1 ML (by ~ 0.56 eV, at large film thickness). This is consistent with the experi- 
mental trend, and an increased CO-CO repulsion at coverages larger than 0.5 MLi2£3i The 
chemisorption energy for the Cu(001)-c(2 x 2) surface (~ 0.74 eV) is consistent with previous 
GGA values,— and somewhat larger than the experimental value (0.57 eV). 2 ^ 



C. Discussion and microscopic interpretation 

In order to better understand the mechanism and identify the key parameter responsible 
for the oscillations in the chemisorption energy, we have investigated the behavior with film 
thickness of (i) the DOS per atom of the films at the Fermi energy, (ii) the work function of 
the films, and (iii) the decay length in vacuum (A) of the electronic local density of states of 
the films at Ep. Both the DOS and A have been invoked as possible key factor responsible 
for the oscillations in the oxidation rate of ultrathin metal films, while the work function is 
a parameter, known to exhibit quantum-size effects,— 1 ^ which has also been been shown to 
be important when discussing reactivity, e.g., in relation with catalytic promotion.— 

In Fig. we display the calculated DOS per atom at Ep of the free-standing Cu(001) 
films, as a function of film thickness. The DOS clearly displays short-period (2 to 3 ML) 
oscillations, which neither follow the periodic crossing between the QWS at T and Ep (in 
Fig. [3]) nor show a correlation with the oscillations in the chemisorption energy of the 
supported films (in Fig. H|). We note that the results in Fig. [5] are consistent with the 
calculated DOS(-Ep) behavior reported for free-standing films with thicknesses larger than 
7 ML in Ref.t. The short-period oscillations of the DOS (Ep) can be related to a beating 
effect due to a supersposition of short and long Fermi-wavelength oscillations.— 

In Fig. El we show the theoretical dependence of the work function on film thickness, for 
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FIG. 5: Calculated density of states per atom at the Fermi energy for the free-standing Cu(OOl) 
films, function of film thickness. 

the free-standing and Fe-supported Cu(OOl) films. The work-function dependence on film 
thickness is similar for the supported and unsupported films, 35 although the amplitude of 
the oscillations is slightly reduced in the supported case. One can notice periodic cusps, 
corresponding to local minima in the work function at 5 ML and 11 ML, which coincide 
with the thicknesses of the crossing of Ep by a QWS at T. These periodic cusps in the work 
function, when a QWS at T crosses Ep, are consistent with the jellium-model predictions by 
Schulte.— We also observe that the work function, in Fig. El is largest at 3 ML and shows 
a local maximum at 8 ML. The calculated value at the largest Cu thickness (4.43 eV at 12 
ML) is close to the experimental work-function value of the Cu(OOl) surface (4.59 eV).— 

An increase/decrease in the work function may be expected to shift the Cu 3d states to 
lower/higher energy with respect to the vacuum level and possibly also with respect to the 
CO electronic levels. Therefore one could expect that maxima/minima in the work function 
may lead to a decreased/increased interaction (hybridization) between the occupied Cu 3d 
states of the surface and the empty CO 2tt* states of the isolated molecule, and hence may 
correspond to minima/maxima in the chemisorption energy— ^ However, comparing the 
oscillations of the work function (Fig. [6]) and of the chemisorption energy of the supported 
films (FigJH), one observes that the local maxima in the chemisorption energy (at 4 and 
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FIG. 6: Calculated work function of the free-standing (left panel) and Fe-supported (righ panel) 
Cu(001) films, as a function of film thickness. 

10 ML) are shifted by -1 ML with respect to the minima in the work function (at 5 and 

11 ML), and that the local minimum of the chemisorption energy (at 6 ML) also does not 
correspond to the local maximum of the work function (at 8 ML). Furthermore, normal- 
emission photoelectron spectral of Cu(001) films indicate a monotonic shift of the Cu 
3<i-levels as a function of film thickness. The 3d centroid moves slightly away from Ep 
with increasing thickness, which would tend to monotonically decrease the Cu(3g?)-CO(27t*) 
interaction and monotonically weaken the Cu-CO bond.— >^ Hence, for CO on Cu films, a 
shift in the 3<i-band energy does not appear to be the dominant factor in determining the 
tend with film thickness of the CO chemisorption energy. 

In Fig. El we present the calculated decay length, A, as a function of Cu film thickness, 
with and without the Fe substrate. The trends are identical in the two cases. The decay 
length exhibits pronounced oscillations, with a first maximum at 5 ML followed by a min- 
imum at 6 ML, a gradual increase to a second maximum at 10 ML, and a subsequently 
gradual decrease for thicknesses up to 12 ML. We note that this behavior exactly coincide 
with the experimental trends of the desorption temperature in Ref. Q. Apart from the local 
maximum at 4 ML, instead of 5 ML, the behavior of the calculated chemisorption energy of 
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the supported films, in FigJU also follows the trends of A in the range 5-12 ML (with a local 
minimum at 6 ML and a local maximum at 10 ML). In Fig. [7J we also reported (see inset) 
the calculated CO chemisorption energy obtained for 1 ML (upper panel) and 0.5 ML (lower 
panel) of CO on the Fe-supported Cu(001) films, as a function of the decay length A of the 
supported films. Except for the cases of 3 and 4 ML — which are the films for which the 
chemisorption energy critically depends on the substrate (see Fig. HJ), there is a rather clear 
correlation (linear relation) between the calculated chemisorption energy and decay length. 
One may note that the variation in the chemisorption energy is rather small (12 meV) in 
the range 5-12 ML. However, this is consistent with the small variation in the desorption 
temperature reported in Ref. [l. 

The oscillations in A are found to be virtually in phase with the periodic crossing of Ep by 
a QWS at f. This is consistent with the predictions for A of a particle-in-a-box model.— ^ 
The model predicts that all states ip n ,k u of a given subband n have the same decay length 



A ra ~ l/y/—E n , where E n is the energy of the subband state n at k\\ = 0, measured relative 
to the vacuum level.— ^ Hence A is dominated by the decay length of the states at Ep 
belonging to the subband whose energy at f is closest to the Fermi level.— With increasing 
width L of the film, the energy of the highest-occupied QWS at f (of the inverted quantum 



well) increases with respect to Ep. Therefore, A first increases as A n ~ l/y/—E n} until the 



QWS n at Y crosses E F , at which point A decreases to the next value A n+ i ~ l/y/—E n+ i 
of the highest-occupied QWS at f ; A then increases again with increasing L, displaying 
systematic oscillations with film thickness L, as observed in Fig. [7J 

The effect of A on the chemisorption energy may be understood in terms of an in- 
creased/decreased overlap, with increased/decreased A, between the most extended QWS's 
of the film (on the vacuum side) at Ep and the frontier orbitals of the molecule,— namely the 
CO 5a highest-occupied molecular orbitals (HOMO) and 2n* lowest-unoccupied molecular 
orbitals (LUMO). The enhanced/reduced hybridization (a) between the unoccupied QWS's 
just above Ep and the CO 5a HOMO and (b) between the occupied QWS's of the film at 
or just below Ep and the CO 2n* LUMO are both expected to strengthen/ weaken the CO 
bonding to the Cu surface.— The former interaction is associated with an electronic charge 
transfer from the 5a orbitals of the molecule to the surface (donation), whereas the latter 
leads to an electronic charge transfer from the metal surface to the CO 2tt* orbitals (back 
donation). The latter interaction is expected to be the dominant one.— 
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FIG. 7: Calculated decay length in vacuum A of the electronic local density of states at the Fermi 
energy for the free-standing (left panel) and Fe-supported (righ panel) Cu(OOl) films, as a function 
of film thickness. The inset shows the calculated CO chemisorption energy obtained for 1 ML 
(upper panel) and 0.5 ML (lower panel) of CO on the Fe-supported Cu(001) films as a function of 
the decay length A of the corresponding Fe-supported Cu film. 

Hence, for Cu thicknesses larger than 4 ML, the trends we obtain in the chemisorption 
energy can be understood in terms of systematic changes in the decay length of the QWS's 
at Ep, which influence the interaction between filled/empty QWS's of the film at Ep~^ •* 
and the LUMO/HOMO of the molecule. We note that the calculated and experimental 
chemisorption energy display the same trend as a function of film thickness in the range 5- 
12 ML. The main differences occurs at 4 ML: while the experimental chemisorption energy 
decreases with decreasing thickness from 5 to 3 ML, the calculated value has a maximum 



14 



at 4 ML. Such a difference, however, can be expected for the thinnest films considered, 
namely 3 and 4 ML, as the behavior of their chemisorption energy sensitively depends on 
the substrate (see Fig. 4). One may therefore also expect the details of the experimental 
interface (e.g., atomic intermixing at the interface) to affect the chemisorption energy of 
these films. 

We note that in our calculations the maximum in the chemisorption energy does not 
occur exactly when a QWS at T crosses Ep, but 1 ML before. We believe this is related 
to the fact that, in the calculations, the quantum- well subband with energy closest to Ep 
at r (see Fig. [2]) barely touches the Fermi level at 5 ML (and 11 ML). In this situation, 
although the unoccupied QWS's just above Ep may interact with the 5a of the molecule, the 
occupied QWS at F cannot interact, by symmetry, with the CO 2n* states. In fact, based 
on projections on the atomic states, we find that the QWS's near Ep at f are composed 
mainly of Cu 4s states, with a small component of Cu 3d a states (about 10 % for the f 
QWS's within 0.2 eV of the Fermi level) M In the region near the Brillouin-zone center, a 
3d„- component, which can mix with the CO(27r*) orbitals, appears in the QWS's, near Ep, 
when going outside f, e.g., along the T — M direction. The 3^ component increases linearly 
with k\\ (reaching, e.g., ~ 10 % at 1/5 of the f - M distance for the 10 ML QWS with 
subband energy at k\\ = closest to Ep in Figj2]). 4 - In order thus for the QWS's with the 
largest possible decay length at Ep to mix with the 2n* states, a small, but non-vanishing 
fraction of the corresponding quantum-well subband should be occupied. This corresponds 
to an optimal situation in which the highest occupied QWS at f is closest to, but not yet 
crossing Ep. In our calculations, this situation corresponds to the 4- and 10- ML cases, 
which in fact yield maxima in the calculated chemisorption energy. 

IV. SUMMARY AND CONCLUSIONS 

We have performed first-principles calculations, based on DFT, to study the modulations 
with film thickness in the chemisorption energy of CO molecules on Cu(001) films. We have 
examined Fe-supported and free-standing Cu(001) films with thicknesses in the range 3-12 
ML. The presence of the substrate has an important influence on the chemsiorption energy 
for thicknesses of 4 ML and below. The calculated CO chemisorption energy of the supported 
films displays systematic oscillations, as a function of film thickness, with a periodicity cor- 
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responding to that of quantum-well states at T crossing the Fermi energy. These oscillations 
in the chemisorption energy are understood in terms of periodic modulations of the decay 
length in vacuum of the quantum-well states at the Fermi energy. These modulations of the 
decay length are expected to influence the interaction of the quantum-well states with the 
frontier orbitals of the molecule. Contrary to previous suggestions, we find that the actual 
oscillations with film thickness of the density of states per atom of the films at the Fermi 
energy cannot account for the observed modulations of the chemisorption properties. 
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